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We present studies of an effective model which is a simple generalization of the standard model 
of a local pair superconductor with on-site pairing (i.e., the model of hard core bosons on a lattice) 
to the case of finite pair binding energy. The tight binding Hamiltonian consists of (i) the effective 
on-site interaction U, (ii) the intersite density-density interactions W between nearest-neighbours, 
and (iii) the intersite charge exchange term I, determining the hopping of electron pairs between 
nearest-neighbour sites. In the analysis of the phase diagrams and thermodynamic properties of this 
model we treat the intersite interactions within the mean-field approximation. Our investigations of 
the U < and W > case show that, depending on the values of interaction parameters, the system 
can exhibit three homogeneous phases: superconducting (SS), charge-ordered (CO) and nonordered 
(NO) as well as the phase separated SS-CO state. 

PACS numbers: 71.10.Fd — Lattice fermion models (Hubbard model, etc.), 74.20.-z — Theories and 
models of superconducting state, 64.75.Gh — Phase separation and segregation in model systems (hard 
spheres, Lennard-Jones, etc.), 71.10.Hf — Non-Fermi-liquid ground states, electron phase diagrams and 
phase transitions in model systems, 71.45.Lr — Charge-density- wave systems 
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I. INTRODUCTION 

The interplay and competition between superconduc- 
tivity and charge orderings is currently under intense 
investigations (among others in high temperature su- 
perconductors such as cuprates, barium bismuthates, 
fullerenes and several other nonconventional supercon- 
ducting materials, e.g., the Chevrel phases) pQ. They 
belong to a unique group of extreme type II supercon- 
ductors and generally exhibit low carrier density, a small 
value of Fermi energy (Ep < 0.1 -j- 0.3 eV) and a short 
coherence length £0 (^o^F ~ 1 10). These general fea- 
tures are consistent with short-range, almost unretarded 
effective interactions response for local pairing. 

In this report we will concentrate on the intriguing 
problem of phase separation and the competition between 
superconductivity and charge orderings. The model 
Hamiltonian considered has the following form: 

H = U^fiifhn - n^fii (1) 
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where ft* = £ CT n«x, h ia = c+cW, pf = (pr)t = g+cj. 
Her (cmt) denotes the creation (annihilation) operator of 
an electron with spin a =t,4- a t the site i, which satisfy 
canonical anticommutation relations: 

{c lCT , ct,} = <%<W,{c iCT , c iff /} = {ct , = 0, (2) 



e-mail: kakonrad@amu.edu. 



where 8^ is the Kronecker delta. Y]u j\ indicates the sum 
over nearest-neighbour sites i and j independently, z will 
denote the number of nearest-neighbours. U, I, and W 
are the interactions parameters, Iq = zl, Wo = zW. p is 
the chemical potential, connected with the concentration 
of electrons by the formula: n = (1/N) J^. {hi), with 
< n < 2 and N is the total number of lattice sites. 

The interactions U, I, and W will be treated as the 
effective ones and will be assumed to include all the 
possible contributions and renormalizations like those 
coming from the strong electron-phonon coupling or 
from the coupling between electrons and other elec- 
tronic subsystems in solid or chemical complexes pQ. 
In such a general case, arbitrary values and signs of 
U, I, and W are important to consider. In the ab- 
sence of the external field conjugated with the SS or- 
der parameter (A = (1/N) ^2i{p7)) there is a symme- 
try between / > (s-pairing) and / < (//-pairing, 77S, 
A Q = (1/N) ^ exp (iQ ■ Ri)(Pi), Q is half of the small- 
est reciprocal lattice vector) cases for model ([lj , thus we 
restrict ourselves only to the / > case. The CO param- 
eter is defined as tiq = (1/N) J^. exp (iQ ■ Ri)(hi) . 

We have performed extensive study of the phase dia- 
grams (PDs) of model ([l]) [2]. In this report we inves- 
tigate the properties of model ([TJ for the case of on-site 
attraction (U < 0, local pair or Bose condensation limit 
[SHE]) and intersite repulsion (W > 0). In the analysis 
of the model at T > 0, we have adopted a variational 
approach (VA) which treats the on-site interaction U ex- 
actly and the intersite interactions W and / within the 
mean- field approximation (MFA). The PDs of model Q 
have been investigated until now for the special cases: 
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FIG. 1. knT / Iq vs. p./Io phase diagrams (upper row) and corresponding ksT/Io vs. n diagrams (lower row) for U/Iq = —1.0 
and different values of W/I = 0.5, 1.1, 2.0 (as labelled). Dotted, solid, and dashed lines indicate first-order, second-order, and 
"third-order" boundaries, respectively. B denotes bicritical points. 



W = [3HT2] and / = [T3H2T] only. Some ground state 
results for W ^ and / ^ have been also obtained 

0|221|23]. 

Within the VA the intersite interactions are decoupled 
within the MFA, which allows us to calculate the averages 
n, tiq, A, and Aq. It gives a set of four self-consistent 
equations (for homogeneous phases). The definitions of 
homogeneous phases with the values of order parameters 
are as follows: (i) SS - n Q = 0, A ^ 0, A Q = 0; (ii) CO 
- n Q £ 0, A = 0, A Q = 0; (iii) M - n Q ^ 0, A ^ 0, 
Aq ^ 0; (iv) NO - riQ = 0, A = 0, Aq = 0. It is 
important to find a solution corresponding to the lowest 
energy. 

Phase separation (PS) is a state in which two do- 
mains with different electron concentration: n + and n_ 
exist in the system (coexistence of two homogeneous 
phases). The free energies of the PS states are calcu- 
lated in a standard way, using Maxwell's construction 
(e.g., Refs. [5J HH OH ). In model Q for the range 
of parameters considered in this paper only one PS state 
can occur, which is coexistence of the SS and CO phases. 

In the report, we have used the following conven- 
tion. A second- (first-)order transition is a transition 
between homogeneous phases with a (dis-) continuous 
change of the order parameter at the transition tem- 
perature. A transition between homogeneous phase and 
PS state is symbolically named as a "third-order" transi- 
tion 0I3H7M]. 



II. RESULTS IN THE LOCAL PAIRING LIMIT 

A. Phase diagrams 

One should noticed that PDs obtained are symmetric 
with respect to half-filling (n = 1) because of the particle- 
hole symmetry of Hamiltonian ([TJ, so the PDs will be 
presented only in the range p, = fj, — U/2 — Wo < and 
< n < 1. 

For any U < and fixed W > 0, the PDs are quali- 
tatively similar, all (first-order, second-order and "third- 
order") transition temperatures decrease monotonically 
with increasing U and in the VA for U — the transi- 
tion temperatures account for a half of those in the limit 
U — > — oo, what can be symbolically written as 

k B T c (U^-^)^2k B T c (U = 0) 1 (3) 

where T c denotes the transition temperature (which can 
be SS-NO, SS-CO, CO-NO, PS-CO or PS-SS). 

Notice that in the U — oo limit model ([!]) is equiv- 
alent with the hard-core boson model on a lattice |25f - 
|2"5] , Moreover in that limit model (|T|) can be derived as 
effective Hamiltonian at the strong-coupling limit of the 
extended Hubbard model by the degenerate perturbation 
theory [231231130]. 

In the range of the attractive on-site interaction the 
structure of PDs of model depends on the ratio W/I 
only (cf. Fig [l] for U/I = -1). One can distinguish 
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two ranges of the ratio W/I in which the system exhibits 
substantially different behaviours: 

(i) < W/I < 1. In Figs. [I^b, we present particu- 
lar PDs for W/I = 0.5. For < W/I < 1 and U < 0, 
only the second-order SS-NO transitions occur with in- 
creasing temperature. If we analyze the system for fixed 
n, for this range of model parameters, the PS states do 
not occur and the obtained PDs have the same struc- 
ture as those derived in [3H5]. The transition between 
homogeneous SS and NO phases taking place with in- 
creasing temperature is second order for arbitrary ft and 
n, and it decreases monotonically with increasing |/2|/io 
and |1 — n\. 

(ii) 1 < W/I. A few particular PDs in this regime are 
presented in Figs.[IJ>f. For W/I = 1 the SS, the CO and 
the M phases are degenerate at n = 1. For W/I > 1, 
three homogeneous phases (SS, CO, NO) appear. The 
SS-NO and CO-NO transitions are of the second order 
and these transition temperatures are decreasing function 
of \p,\ and |1 — n\. The SS-CO transition is discontinuous 
for fixed p,, and thus the PS state SS-CO is stable in the 
definite range of n. All transitions lines meet at a bicriti- 
cal point B. With increasing W, the B-point moves along 
the boundary between SS and NO phases toward larger 
|/2 1 (|1 — n|). This is due to the fact that in the VA the 
SS-NO transition is independent of W (for fixed n). The 
region of the CO phase occurrence is extended, whereas 
the region of the SS phase stability is reduced by increas- 
ing the ratio W/I. The first-order SS-CO as well as the 
"third-order" SS-PS and PS-CO transition temperatures 
increase with \p,\ and |1 — n|, respectively. 

One should notice that for W^ restricted to nearest 
neighbours {W% = 0) the PS state is strictly degenerated 
at T = with the M phase in the whole range of sta- 
bility of both these states |22[ [23] . This degeneration is 
removed at T > 0, even for W% = and the PS state 
occurs on PDs. Repulsive W2 > between next-nearest 
neighbours destabilizes the PS state with respect to the 
M phase, whereas attractive W2 < extends the stability 
region of PS state and eliminates the M phase. 

In the case of attractive W < (precisely for W/I < 
— 1), the model can exhibit phase separation NO-NO 
(electron droplets state) at low temperatures |201 l2"2l 12"51 
I3T] and for U < the PDs as a function of n have the 
similar structure as those derived in [!JT] . We leave deeper 
analysis of this problem to future publications. 



B. Order parameters and thermodynamic 
properties 

Let us focus now on the temperature dependencies of 
the order parameters and thermodynamic properties of 
the system at the sequence of transitions: SS — > CO — > 
NO for W/I = 2.0, U/Iq = -1.0 and fi/I = -2.0. 

The temperature dependencies of the order parame- 
ters: A and nq are presented in Fig. [2] It is clearly 
seen that at the SS-CO transition (at fcsT c i//o = 0.40) 
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FIG. 2. Temperature dependencies of (a) superconducting 
order parameter |A| and (b) charge order parameter uq for 
W/I = 2.0, U/I = -1.0 and fl/I = -2.0. 



the both order parameters change discontinuously. The 
CO-NO transition at ksT^/Io = 0.83 is of the second 
order. 

Calculating the grand potential per site uj = 
-l/(iV/3)ln(Tr[exp(-/3.H r )]) within the VA one can ob- 
tain thermodynamic characteristics of the system for ar- 
bitrary temperature. The entropy s and the specific 
heat c can be derived as s = —(duj/dT) and c = 
—T(d 2 uj/dT 2 ). uj, s, and c as a function of temperature 
are shown in Fig. [3] s increases monotonically with in- 
creasing T. It is discontinuous at T c \ whereas it is contin- 
uous at T C 2- One can notice that in the high-temperature 
limit s/ks —> hi(4) w 1.386 (there are four possible con- 
figurations at each site). The peak in c(T) is associated 
with the first-order transition (at T c i), while the A-point 
behaviour is typical for the second-order transition (at 
T c2 ). 



III. CONCLUSIONS AND FINAL REMARKS 

In this report, we have presented some particular PDs 
of the extended Hubbard model with pair hopping and in- 
tersite density-density interactions in the zero-bandwidth 
limit for the case of local attraction U < 0. One finds that 
the system considered can exhibit very interesting mul- 
ticritical behaviours. Our investigations show that, de- 
pending on the values of interaction parameters (the ratio 
W/I), the system can exhibit three homogeneous phases: 
superconducting, charge-ordered and nonordered. The 
SS-NO and CO-NO transitions are of the second order. 
The SS-CO transition is discontinuous (for fixed fx) , what 
leads to phase separation on the phase diagrams as a 
function of n for W/I > 1. The homogeneous mixed 
phase (with nonzero both charge-ordered and supercon- 
ducting order parameters) never occurs on PDs at T > 0, 
at least in the absence of the next-nearest neighbours in- 
teractions. On the contrary the PS state: SS-CO is found 
to be stable in definite ranges of model parameters and 
temperatures. 

Our results are exact in the limit of infinite dimensions. 
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FIG. 3. Temperature dependencies of thermodynamics parameters: (a) the grand potential to, (b) the entropy s, and (c) the 
specific heat c for W/I = 2.0, U/I = -1.0 and p/I = -2.0. 



where the MFA treatment of intersite terms is rigourous 
one [TT] . In finite dimensions due to quantum fluc- 
tuations connected with the I term, the regions of the 
ordered homogeneous phases occurrence are extended in 
comparison with the VA results [B]. 
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